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Abstract 

The objective of this paper is to investigate the stability of limit cycles 
of a mathematical model with a distributed delay which describes the 
interaction between p53 and mdm2. Choosing the delay as a bifurca- 
tion parameter we study the direction and stability of the bifurcating 
periodic solutions using the normal form and the center manifold the- 
orem. Some numerical examples are finally made in order to confirm 
the theoretical results. 
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1. Introduction 

The tumor suppresser gene p53 plays a key role in oncogenesis and its 
anomalies are almost universal in tumoral cells [4]. To make the reading 
easier we briefly present the role of p53 gene in living cells. Normally the 
activity of p53 is inhibited and the concentration of p53 protein are kept at 
very low levels. The activation occurs when there are DNA damages [2,3]. 
Depending on these damages there are two outcomes: one is the cell cycle 
arrest induced by a low level or a brief elevation of p53 protein, and the 
other is the apoptosis induced by a high level or a prolonged elevation of p53 
protein [5]. Knowing these 2 different outcomes, it is now clear that the level 
of p53 protein should be kept tight control. This control is achieved with the 
help of mdm2 gene with which p53 makes a feedback loop [5, 10]. Recently it 
had been discovered that this loop is not quite straightforward because there 
are two isoforms of the mdm2 protein, i.e. p76mdm2 and p90mdm2 which 
have different features and roles [7]. 

This new model is based on [8, 9]. Here we achieve a smoother modeling 
of the phenomenon, i.e. the interaction p53-mdm2. The production of p53 
protein is continuous, so is the binding between p53 and the promoter of the 
mdm2. The difference from the previous model [9] lies in the introduction of 
the integral form in the third equation, which is the natural way of modeling 
a continuous process. To have a clear picture of this aspect, we explain the 
need for the introduction of the integral form, by the fact that the synthesis of 
p53 protein and the binding p53-mdm2 are continuous and there is no buffer 
to store p53 and then, at a given moment, to release it, and after that all the 
quantity of p53 to bind mdm2 promoter. As a matter of fact, the molecules 
of p53 protein that are bound with mdm2 promoter were not synthesized at 
the moment t, but they were synthesized at a previous different moments, 
they were bound at different moments. 

The variables of the model are: xi, X2 mRNA concentrations and yi, 2/2 
the protein concentrations. 
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The mathematical system which describes our model is: 



= 1 - hx^t), 
ijiit) = - (ai + a 12 y 2 (t))yi(t), 

1 r (!) 

x 2 (i) = - / (af( yi {t)) + (1 - a / j/i t - s)))ds - b 2 x 2 (t), 

T JO 

y 2 (t) = x 2 (t) - (a 2 + a 12 yi{t))y 2 (t) 

where: bi, b 2 are the rates for mRNA degradation, a±,a 2 , a\ 2 are the rates for 
protein degradation. The function / : 1R + — > 1R, is the Hill function, given 
by: 

f(*) = ^T—n (2) 

a + x n 

with n G IN*, a > 0. The parameters a±, a 2 , b±, b 2 , a i2 of the model are 
assumed to be positive numbers less or equal to 1, a G [0, 1] and r > 0. 

For a = in the present model, we obtain the model from [9], which 
suggests that there is an oscillatory behavior based on using only numerical 
simulations. 

For the study of the model (1) we consider the following initial values: 

xi(0) = xi,yi(0) = (pi(6),6 G [-t,0],x 2 (0) = x 2 ,y 2 (0) = y 2 , 

with x\ > 0, x 2 > 0, y 2 > 0, (p\(6) > 0, for all 6 G [— t, 0] and ip\ is a 
differentiable function. 

The paper is organized as follows. In section 2, we discuss the local sta- 
bility for the equilibrium state of system (1). We investigate the existence 
of the Hopf bifurcation for system (1) using time delay as the bifurcation 
parameter. In section 3, the direction of Hopf bifurcation is analyzed by the 
normal form theory and the center manifold theorem introduced by Hassard 
[6]. Numerical simulations for justifying the theoretical results are illustrated 
in section 4. Finally, some conclusions are made and further research direc- 
tions are presented. 

2. Local stability and the existence of the Hopf bifurcation. 

The equilibrium points of system (1) are given by the solutions of the 
following system of equations: 
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1 - = 0, 

xi(t) - (ai + a 12 y 2 (t))y 1 (t) = 0, 

x 2 (t) - (a 2 + a 12 yi(t))y 2 (t) = 0. 
Let g : (0, oo) — > IR be the function given by: 



(3) 



(a 2 + ai 2 x)(l - aiftix) , , 

bidi 2 X 

From (3) and (4) it results that one solution of system (3) is: 

1 (a 2 + a 12 y 1 )(-a 1 b 1 y 1 + 1) -a 1 b 1 y 1 + 1 
xi = 7-, yi, x 2 = , y 2 = — , (5) 

&i &iai2j/l Ml2?/l 

where 2/1 is the solution of the equation: 

f(x)-g(x)=0. (6) 

Because /(x) is an increasing function with lim /(x) = 1 then g(x) is an 

increasing function too, for x < . It results that equation (6) has a 

ai&i 

unique solution. 

Proposition 1. //yio is a real solution of (6) then the equilibrium point 
of system (1) is: 

1 1 - ai&iyio (a 2 + ai 2 yi )(l - ai&i2/io) 

£io = 7-, J/20 = —7 , ^20 = r • (7) 

h Mi2j/io M122/10 

We consider the following translation: 

xi = Mi + xio, 2/i = u 2 + j/io, x 2 = w 3 + x 20 , y 2 = u 4 + 2/20- (8) 
With respect to (8), system (1) can be expressed as: 
ui(t) = -6itii(t), 

«2(*) = - (ai + a 12 y 20 )u 2 (t) - a 12 y 10 u 4 (t) - a 12 u 2 (t)u 4 (t) , 

i r ( 9 ) 

«3(*) = - / («/( M 2(0+2/lo) + (l-a)/(M2(^-s)+l/lo))^-&2M3(t)-&2^20, 

T JO 

u 4 (t) = -a 12 y 20 u 2 (t) +u 3 {t)-(a 2 + a 12 y w )u 4 {t)-a 12 u 2 {t)u 4 (t). 
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System (9) has (0,0,0,0) as equilibrium point. 

To investigate the local stability of the equilibrium state we linearize 
system (9). We expand it in a Taylor series around the origin and neglect 
the terms of higher order than the first order for the functions from the right 
side of (9). We obtain: 



where 



U{t) = AU(t) + -B / U(t - s)ds, 



(10) 



A= 



f-h 
1 
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-{a 2 +a 12 y w ) 
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(11) 



with Pl = f(y w ), u(t) = { Ul {t),u 2 {t)Mt)Mt)) T , foU(t-s)ds = UoMt- 

s)ds,f G T u 2 (t - s)dsJv u 3 (t - s)ds, J r u 4 (t - s)ds) T . 

The characteristic equation corresponding to system (10) is det(XI — A — 



T JO 



e~ Xs ds)B) = which leads to: 



(A + foiXA 3 + b\ 2 + cA + d + - f T e- Xs ds) = 0, 

r Jo 



(12) 



where 



b = oi + a 2 + b 2 + a 12 (y 20 + y 10 ), 
c = b 2 {ai + a 2 ) + b 2 a 12 (y 20 + y w ) 
d = b 2 a x a 2 + a 12 b 2 (y 20 a 2 + aiy 10 ) + aa 12 y 10 p 1 , 
h = (1 - a)a 12 y w pi. 



a x a 2 + ai 2 (aiyi + 02^/20), 



(13) 



The equilibrium point X* = (a: 10 , yio, ^20? U2o) T is locally asymptotically 



stable if and only if all eigenvalues of (12) have negative real parts. 
Because b\ > we will analyze the function: 



A(A, r) = A 3 + b\ 2 + cA + d + - / 

T JO 



-As 



ds, 



(14) 
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with A G H. 

We are going to show that the equilibrium point X* undergoes a Hopf 
bifurcation. In this sense, we look for the existence of the purely imaginary 
roots of A(A, r) = 0. First, we verify if X* is locally asymptotically stable 
when r = 0. In this case, the equation A(A, r) = becomes: 

A 3 + b\ 2 + c\ + d + h = 0. (15) 

Because the coefficients of equation (15) are positive then according to the 
Routh-Hurwitz criterion we have: 

Proposition 2. When there is no delay, the equilibrium point X* of 
system (1) is locally asymptotically stable if and only if 

cb > d + h, 

where c, b, d, h are given by (13). 

We are looking for the values To so that the equilibrium point X* changes 
from local asymptotic stability to instability or vice versa. This is specific for 
the characteristic equation with pure imaginary solutions. Let A = ±iu be 
these solutions. We assume uj > 0. It is sufficient to look for A = iu root of 
A(A, r) = 0. Separating real and imaginary parts of A(iu, r) = we obtain: 

. , , Tuj(buj 2 - d) . , _ tuj 2 (c-uj 2 ) , . 

siu{ujt) = — y — cos{ujt) = 1 y — -. (16) 

lb lb 

A solution of (16) is a pair (uj,t) so that cut G [0,27r], uj G (0,y/c\ and 
r = gi(co), where g± : [0, \/c\ — ► 1R + is given by: 

9l{X) ~ (bx 2 -df + x 2 (c-x 2 f { } 

2hc 2hc 

From (17) it results that gAx) G [0, —pA and for all r G [0, — ^-] there is 

d 2 d 2 

uj G [0, ^/c\ so that r = g(u>). 
We have: 

2hc 

Proposition 3. For t G [0, -^-] there is uj G [0, -^/c], so that (uj,t) is 
solution of (16). 

According to the above proposition, system (16) has not unique solution. 
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In order to show that X* undergoes a Hopf bifurcation for r = tq wc 
have to prove that ±iuo are simple eigenvalues of A (A, To) and satisfy the 

dRe(\). 

transversahty condition T=T „ 7^ 0. 

dr 

From (14), it results that: 

A A (A, r) = |±(A, r) = 3A 2 + 2b\ + c - -£-(1 - (1 + Ar)e^)- (18) 

From (18), for (cu,r) a solution of (16) and r = g(u;) it results that: 

M l = Re(A\(iu, gi(u)) = -Au 2 + 2c- r{buj 2 - d) 

T /a /■ /xx 3buj 2 -d-h 2 , 

M 2 = lm(A\(iuj } gi{u!)) = h ro;(c -w J. 

(J 

Then, we have: 

M 2 +M 2 = (2c-4u 2 ) 2 + T 2 (bLU 2 -d) 2 + {3buj2 ~^ h)2 +t 2 u 2 (c-uj 2 ) 2 + 
+ 2r({3buj 2 -d-h)(c-u 2 ) -{bu 2 -d){2c-4uj 2 )). 

(19) 

From (19), it results that: 

Proposition 4. If be > d + h, d > h, be > 3d — h, then A = iujq is a 
simple root for the equation A(A, To) = 0. 

Now, we consider a branch of eigenvalues A(r) = v{t) + icu(r) of (14) so 
that v(to) = and uj{tq) = u , where (cj ,to) is a solution of (16). Differen- 
tiating equation A(A(r),r) = with respect ro r, we obtain: 

\'( \ dX \ A T (iw ,r ) \ . • at/ v 

A (r) = —\ T=T0!U)=U , = - 1 — r . r = M{uq, t ) + iN(v , r ) 

dr A A (^o,r ) 



where 



with 



and 



M 1 N 1 + M 2 N 2 
M( "°' ro) = - M 2 + Mf (20) 

iVx = -[/i - (k^ 2 - d) - tu 2 (c - u 2 )} 

I 

N 2 = -[uj{c - u 2 ) - Tuj{buj 2 - d)] 

T 

AT M,N 2 - M 2 N, 

7V " M 2 + M 2 ■ [Zi) 
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By direct calculation M^O. 

2hc 

We can conclude that when r G [0, — — ] the characteristic equation 

cr 

A(A, ro) = has a unique pair of purely imaginary simple eigenvalues satis- 
fied Q ^ e ^( r ) | j. q an j foe real roots are negative. Consequently a Hopf 
(It 

bifurcation occurs at X* when t = tq. Moreover, applying Rouche's theorem 

we can verify that every eigenvalue of A(A, r) = with r < tq has negative 

real part. It follows that X* is locally asymptotically stable for < r < r . 

These results are summed up in the following theorem: 

Theorem 1. Assume that cb > d + h. Then there exists values r G 
2hc 

[0, — ^-] o/ time c?eZay so £/ia£ £/ie equilibrium point X* is locally asymptoti- 
cally stable when r G [0, r ) and becomes unstable when r = r throughout a 
Hopf bifurcation. In particular, the periodic solutions appear for system (1) 
when t = t . 



3. Direction and stability of the Hopf bifurcation 

In the previous section, we obtain some conditions with guarantee that 
system (1) undergoes Hopf bifurcation at r = r . 

In this section, we study the direction, the stability and the period of the 
bifurcating periodic solutions. The used method is based on the normal form 
theory and the center manifold theorem introduced by Hassard [6]. 

For an interval / C ]R, we define the space of continuous functions as 
C(I,K) = {/ : / -> K,f continuous}, where K = H 4 or(D 4 . When / = 
[— r, 0], r = tq + fx, \x > sufficiently small, we set C M = C([—r,0],K). 
Expanding the functions from the right side of system (9) in Taylor series 
around (0, 0, 0, 0) T it results that: 

X(t)=AX(t)+-B f X(t- s)ds+F(X(t), f \aX(tyAl-a)X(t-s))ds) (22) 
r Jo jo 

where 

X(t) = ( Ul (t), U2 (t),u 3 (t), U4 (t)) T , 

f'T PT f'T PT PT 

/ X(t—s)ds — ( u 1 (t—s)ds, U2(t—s)ds, us(t—s)ds, U4(t—s)ds) T , 
jo jo jo jo jo 



8 



F(X(t) ,£(aX(t) + (1 - a)X(t - s))ds) = (0,F 2 (« 2 (t),« 4 (t)), 
F\^{au 2 {t) + (1 - a)u 2 (t - s))ds),F\u 2 (t), u 4 (t))) T , ( 23 ) 

F 2 (u 2 {t), u 4 (t))=-a 12 u 2 (t)u 4 (t), 

F 3 {£{au 2 {t) + (l-a)u 2 {t-s))ds)~p 2 f \au 2 {t) + {l-a)u 2 {t-s)) 2 d^ 

It Jo 

! FT 

+7T-P3 / (au 2 (t) + (1 - a)it 2 (t - s)) 3 ds, 
or jo 

F 4 (w 2 (t),w 4 (t)) = -ai2U 2 (t)u 4 (t), 

P2 = /"(yio), P3 = /'"(2/io) and A,B are given by (11). 
For $ G C M with K =(E 4 we define a linear operator: 



L M ($) = A$(0) - — 5 /° $(s)ds 

Tn J-rii 



To J -to 

and a nonlinear operator: 



F„($) = (0,F 2 ($ 2 (0),$4(0)),F 3 (|% 2 ( S )rf S ),F 4 ($ 2 (0),$ 4 (0))) T . 

For $ G C^Q-rcO]^ 4 ) we define: 

d$(9) 



6e[-T ,0) 



d9 ' 

A$(0) - —5 /° $(s)ds, = 0, 

Tn J-rn 



0, 0G[-r o ,O) 

e = o 

and for \I/ G C 1 ([0, r ],(E* 4 ), we define the adjoint operator A* of A by: 

d*(s) 



„4**(s) 



ds 

l , r° 



s G (0, r ] 



A*(0) + -( P V(O)d0)B, s = 0. 

Tn •'0 



T 

Then, we can rewrite (22) in the following vector form: 

X t = A(n)X t + R(n)X t (24) 
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where X t = X(t + 9) for 6 G [-r ,0]. For $ G C([-r , 0],(E* 4 ) and # G 
C([0, tq],(D* ) we define the following bilinear form: 



< *(s)mo) >= *(o)*(o) - /° / e - /* $(CK'R^ 

J -to J£=0 T JO 



se [0, To], # G [—To, 0]. 

Then, it can verified that ^4* and ^4(0) are adjoint operators with respect 
to this bilinear form. 

In the light of the obtained results in the last section, we assume that 
±iui are eigenvalues of .4.(0). Thus, they are also eigenvalues of A*. We can 
easily obtain: 

§(9)=ve Xl \ 0e [-7b,0] (25) 
where v— (vi,v 2 ,v 3 , v 4 ) T , 

vi = 0,v 2 = a 12 y w , v 3 = a 2 12 y w y 20 - (X 1 + en + a 12 y 20 )(X 1 + a 2 + a 12 y 10 ), 
v A = -(Ai + ai + a 12 y 20 ) 

is the eigenvector of .4(0) corresponding to Ai = iu> and 

q(s)=we XlS , sG[0,r ] 

where w = (wi, w 2 ,w 3 , w±), 

1 d 2 d 3 di 
w 1 = -,w 2 = —,w 3 = —,w 4 = —, 
7] 77 rj 77 

r , , , , ai2j/io(6i + Ai) ai 2 yio(6i + Ai) 

d 2 = Oi + Ai, d 3 = -y- — — , d 4 = 



(Ai + a 2 + ai 2 yio)(o 2 + Ai) Ai + a 2 + a 12 y 10 

7] = v 2 d 2 + v 3 d 3 + vj 4 - ( W 1 " 3 1 (r Ai - 2 + 2e- Airo + AiTbe"* 176 ) 

To Ax 

is the eigenvector of .4* corresponding to A 2 = — iu> - 

We can verify that: < #(s),$(s) >= 1, < *(s),l>(s) >=< *(s), <&(s) >= 
0, < *(s),l>(s) >= 1. 

Using the approach in [1], we next compute the coordinates to describe 
the center manifold f2 at ii = 0. Let X t = X(t + 9), 9 G [— r ,0], be the 
solution of equation (24) when ii — and 

z(t) =< AT t >, u;(t, 0) = X t (0) - 2 J Re{z(t)$(6>)}. 
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On the center manifold f2 , we have: 

w(t,9) = w(z(t),z(t),9) 

where 

w(z, z, 0) = w 20 (9) — + w u (9)zz + w 02 (9) — + w 3Q (9) — + ... 

in which z and z are local coordinates for the center manifold Qo i n the 
direction of \I/ and \l/ and wo2(9) = W2o(9). 

For solution X t of equation (24), as long as fj, — 0, we have: 

i(t) = Ai2(t)+ 5 (z,z) (26) 

where 



z) = V(0)F(w(z, z, 0) + f?e(z$(0))) = 

z 2 z^z 
= 920— + guzz + g 02 — + 92i— + • • • 

From (23), (26) and (27) we obtain: 

#20 = ifo^ + ^lo^3 + F£ u>4, 9u = F^w 2 + F^w 3 + F^w 4 , 
902 = F* 2 w 2 + F$ 2 w 3 + F 4 2 w 4 , 

where 



(27) 



(28) 



F 2D = F£ = -2a l2 v 2 v A , = = -a 12 {v 2 v 4 + v 2 v A ), 

72 - Z?4 
02 — r 02 



F m = FL = -2a 12 v 2 v 4 



Flo = ^f(2a 2 r A 1 - 4a(l - a){e~ x ™ - 1) - (1 - af{e~ 2X ^ - 1)) 

F 3 = P2^2 (AiA2(ft2 + (1 _ a)2)ro _ _ a)(A2e -A lT0 + Aie -A 2 r 0)) 

r AiA 2 

= ^l(2a 2 T \ 2 - 4a(l - a)(e^ - 1) - (1 - a) 2 (e" 2A2 - - 1)) 



and 



92i = F*w 2 + F*w 3 + F*w A (29) 
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where 



F ii = F 2i = -ai2^2«4o( ) - 2a 12 v 2 wf 1 (0) - 012^4^20(0) - 2a 12 v A wl 1 (0) 
F 2 \ = ^(2v 2 (a\wU0) + a(l - a)h - "^f^^ e^ - 1) + 

'Mi Ai 

+ (1 - afk 2 ) + 2v 2 (a 2 r w 2 20 (0) - ° (1 ~ a) <(0)(e" AlT ° - 1) + 

A2 

+a(l - a)fc 3 + (1 - a) 2 A; 4 )) + 

+ ^^^( a 3 To _(l^ (e - 2 A 1 ._ l) _(l^ (e -A 2 ._ lWl ^ ) 3 To) 
r 2Ai A 2 

with 

h— w 2 1 (-s)ds,k 2 = / e~ Xls w 2 1 (-s)ds, 
Jo Jo 

ks = wl (—s)ds,ki= / e~ X2S wl (—s)ds 
Jo Jo 

wU-s) = -fv 2 e-^ - ff^ 2 e- A - + E\e^ 

wK-s) = 9 -^v 2 e'^ - ^vle-^ + E\ 
\\ Ai 

2 / n \ ^20 #02 - . j-,2 

^20(0) = -— v 2 - —v 2 + E 2 

2 m\ 9n #11 - . 
w ii (0) = -T-V2 - -T-V2 + E l 
Ai A]_ 

4 / n x #20 #20 - . T-i4 

^20(°) = ~~ U4 - ^X V ± + E 2 

4 / n \ # n #11 - . 
<i(°) = ^-^4 - -T- v 4 + E u 

Ai Ai 

s E [0, To], E\ respectively E\, E 2 are the components of the vectors: 

E 2 = -(A - J—( e - 2X ^ - l)B - 2A 1 /)- 1 F 20 
2AiT 

E 1 = -(A + B)- 1 F 11 , 

where F 20 = (0, F 2 , F 2 3 , F 2 4 ) T , F u = (0, F*, F*, F^f. 

Based on the above analysis and calculation, we can see that each in 
(28), (29) is determined by the parameters and delay from system (1). Thus, 



(30) 
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we can explicitly compute the following quantities: 



Ci(0) = 2^(^20511 - 2|£ii| 2 - ^o 2 | 2 ) + 9 -f 

Re?CM) 1171(0,(0)) + ^N(u ,r ) 

112 = = ^ ' = 2i?e(Cl(0)) ' 

(31) 

where M(u ,t ), N(uq,t ) are given by (20) and (21). 
In summary, this leads to the following result: 

Proposition 5. In formulas (31), fj, 2 determines the direction of the 
Hopf bifurcation: if /i 2 > 0(< 0), then the Hopf bifurcation is supercritical 
(subcritical) and the bifurcating periodic solutions exist for r > r (< r ); (3 2 
determines the stability of the bifurcating periodic solutions: the solutions are 
orbitally stable (unstable) if [3 2 < 0(> 0); andT 2 determines the period of the 
bifurcating periodic solutions: the period increases (decreases) ifT 2 > 0(< 0). 



4. Numerical examples. 

For the numerical simulations we use Maple 9.5. In this section, we 
consider system (1) with a 1 = a 2 = 0.13, a 12 = 0.02, a 12 = 0.06, &i = 0.2, 
b 2 = 0.4, a = 4, a = 0.2, n — 3. Waveform plot are obtained by the formula: 

X(t+9) = z(t)^9)+z(t)^6}^w 2 o(9)z 2 (t)+w 11 (d)z(t)^ 

where z(t) is the solution of (26), $(6>) is given by (25), w 20 (6), Wn(9), w 02 (9) 
are given by (30) and X = (x±o, yw, x 2 o, y 2 o) T is the equilibrium state. 

We obtain: x 10 = 5, y w = 21.03417191, y 20 = 1.795140515, x 20 = 
2.498925919, fi 2 = -0.2101567953, (3 2 = -0.3029980114, T 2 = 0.1148699183, 
a; = 0.1, r = 0.1001651263. Then the Hopf bifurcation is subcritical, the so- 
lutions are orbitally stable and the period of the solution is increasing. The 
wave plots are displayed in figl and fig2: 
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Fig.l. (t, yi (t)) 


Fig.2. (t,y 2 (t)) 
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5. Conclusions. 

As in our previous model [9], we obtain an oscillatory behavior similar 
to that observed experimentally [3]. The conclusion is not surprising, but 
is useful as this model provides a more accurate approach of the interaction 
p53-mdm2. We can conclude that the transformation made by us to the 
continuous model with distributed time of the interaction p53-mdm2, which 
actually is more real, did not alter the behavior of the system. 
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